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We study discrete surface solitons in semi-infinite, one-dimensional, nonlinear (Kerr), quasiperi- 
odic waveguide arrays of the Fibonacci and Aubry-Andre types, and explore different families of 
localized surface modes, as a function of optical power content ('nonlinearity') and quasiperiodic 
strength ('disorder'). We find a strong asymmetry in the power content of the mode as a function 
of the propagation constant, between the cases of focussing and defocussing nonlinearity, in both 
models. We also examine the dynamical evolution of a completely-localized initial excitation at the 
array surface. We find that in general, for a given optical power, a smaller quasiperiodic strength is 
required to effect localization at the surface than in the bulk. Also, for fixed quasiperiodic strength, 
a smaller optical power is needed to localize the excitation at the edge than inside the bulk. 
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I. INTRODUCTION 

Optical waveguide arrays and photonic lattices have 
proven particularly useful to the study and observation 
of the general phenomenology of wave localization aris- 
ing from the interplay of nonlinearity, periodicity (or 
lack of it), disorder, boundary effects, geometry and di- 
mensionality. Among other phenomena, it has permit- 
ted the direct observation of Anderson localization [T], 
Tamm states in truncated optical lattices [2] and bound- 
ary effects on Anderson localization j3]. Regarding the 
effects of quasiperiodicity in photonic systems, low-loss 
resonant modes have been predicted in aperiodic nanopil- 
lar waveguides [3] , and localization of interface modes at 
the boundary between metal and Fibonacci quasiperiodic 
structure have been recently studied [5]. Also, for a one- 
dimensional quasiperiodic waveguide array known as the 
Aubry-Andre (AA) model, a localization transition was 
experimentally observed |6] in agreement with the ear- 
lier prediction [7 . Moreover, and in contrast to periodic 
systems, quasiperiodic lattices support the existence of 
linear localized surface modes [8] and a rich hierarchi- 
cal band-gap structure [SJ. For this reason the region of 
existence (in parameter space) of surface gap solitons is 
expected to be larger in the quasiperiodic case than in 
the periodic case. 

Concerning the interplay of boundary and disorder, 
Szameit et. al.p] observed experimentally that, in a 
truncated disordered waveguide array, localization due 
to disorder is weakened near the boundary of the trun- 
cated array. This result, in addition to the well-known 
fact that nonlinear localization at the edge of a trun- 
cated periodic nonlinear (Kerr) lattice is also weakened 
by the presence of a boundary [TU], reinforces the idea 
that in one dimension, the boundaries are "repulsive" . 
The whole story is not that simple, however. A recent 
work[TT] extended Szameit's studies by exploring the in- 
terplay of boundaries, disorder and nonlinearity, finding 
that the boundary is "repulsive" or not depending on 
the relative strength of nonlinearity and disorder: For 
weak nonlinearity and moderate disorder, localization is 



weaker at the surface than in the bulk, in agreement with 
Szameit's results. However, for relatively strong disorder 
and/or nonlinearity it is as easy to localize an excitation 
at the edge as in the bulk. 

In this work we explore another side of these issues by 
replacing in these studies the presence of uncorrelated 
disorder by the much weaker loss of periodicity, produced 
by quasiperiodicity. That is, we examine the interplay 
of boundary effects, quasiperiodicity and nonlinearity, in 
two different semi-infinite nonlinear optical waveguide ar- 
rays. In one of them, the linear index of refraction follows 
the Fibonacci sequence, while in the other one, the in- 
dex of each fiber varies in space in an incommensurate 
manner, according to the AA model. The strength of 
quasiperiodicity is varied by adjusting the relative re- 
fractive index contrast of the guides, while nonlinearity 
is regulated by the power content of the input beam. 

We find different families of discrete solitons with dif- 
ferent power thresholds: High power solitons associated 
to total internal reflection (TIR) gaps with finite power 
threshold, and low power solitons associated to Bragg re- 
flection (BR) gaps where the threshold value depend on 
the symmetry of linear modes and their propagation con- 
stants. Contrary to what happens in periodic systems, 
here we have localized linear modes and extended asym- 
metric linear modes that can originate discrete solitons 
with zero power threshold. For a fixed distribution of 
refractive indices, we find a strong asymmetry between 
the focussing and defocussing cases, when the quasiperi- 
odicity of the array is diagonal, i.e., when it affects the 
index of refraction of the guides but not their mutual 
coupling. We also examine the dynamical evolution of 
a completely-localized initial excitation at the array sur- 
face. We find that in general, for a given power content, a 
smaller quasiperiodic strength is required to effect local- 
ization at the surface than in the bulk. Also, for a fixed 
quasiperiodic strength, a smaller power is needed to lo- 
calize the excitation at the edge than inside the bulk. 
This behavior is in marked contrast to the case of uncor- 
related disorder and weakens the usual assumption that 
a surface acts in a 'repulsive' manner [51 110] . 
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This paper is organized as follows. In Sec. II we de- 
scribe the two kinds of quasiperiodicity (substitutional 
and incommensurate arrays) and we discuss the appro- 
priate adimensionalization in each case. In Sec. Ill we 
study the different families of surface localized modes and 
their linear stability. In Section IV we examine the inter- 
play between quasiperiodicity, boundary and nonlinear 
effects on the dynamical evolution of an initially localized 
input beam and we compared the results obtained for sur- 
face and bulk excitation. In section V we discuss briefly 
the nondiagonal quasiperiodic case and finally, section VI 
concludes the paper. 



II. MODEL 

Let us consider a semi-infinite, weakly-coupled, one- 
dimensional quasiperiodic array of nonidcntical nonlin- 
ear optical single-mode waveguides. In the coupled-mode 
framework, the electric field E(x) propagating along the 
array can be presented as a superposition of the modes 
of each waveguide, E(x) = ^ n E n <j) n {x) , where E n is 
the amplitude of the (single) guide mode <j) n {x) centered 
on the nth site. The evolution equations for the modal 
amplitudes E n take the form 



dz 



e n E n + V{E n+1 + E n _{) + 1 \E n \ 1 E n =0 (1) 



where n denotes the position of a guide center, z is the 
longitudinal propagation coordinate, 7 is the effective 
nonlinear coefficient, V is the coupling between nearest- 
neighbor guides and e n is the refractive index of the n- 
th guide. Equation ([lj is known as the Discrete Non- 
linear Schrodinger (DNLS) equation and is also found 
in the study of intrinsic localized modes [12] , molecular 
crystals [15] . biopolymers [13] . arrays of Josephson junc- 
tions [15], Bose- Einstein condensates in magneto-optical 
traps [IB] , among others. 

In this work the sequence {e n } will not be random but 
aperiodic. We will consider two kinds of quasiperiodic 
arrays: substitutional (Fibonacci) and incommensurate 
(AA). 



A. Substitutional arrays 
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FIG. 1. (Color online). Linear (7 = 0) spectrum as a function 
of e for Fibonacci (top) and Aubry-Andre (bottom) arrays. 
For the AA model, x = (v5 + l)/2. (a) Fibonacci surface 
mode, (b) and (c) show an asymmetric extended mode for 
the AA model, before (e = 1) and after (e = 3) the local- 
ization transition (at e = 2), respectively, and for the same 
eigenvalue. 



nonlinearity (7 = 0), the linear modes are described by 
A n (z) = A n e lXz , reducing Eq.|l]) to an eigenvalues prob- 
lem. In the limit N — ¥ 00, the linear spectrum converges 
to a fractal structure, with two main bands (Fig.l), each 
of which is in turn endowed with an internal, self-similar 
structure of bands and gaps [17] . The eigenstates has a 
self-similar structure too in accordance with the aperi- 
odicity of the array. Among these states, there are an 
infinite number of surface localized modes. For a finite 
array, and depending on the value of N, the number of 
such linear surface states is finite, appearing as isolated 
levels of the linear spectrum. It is not easy to determine 
how many of these linear surface states exist for a given 
N. For instance, for an array with N = 89 waveguides, 
we find at least one surface localized mode (Figs, la, 2b). 



The first case corresponds to a substitutional array 
formed by the successive application of the Fibonacci 
substitutional rule: {e a — ¥ e a £&, £b £a}, starting from 
€q = e a . This leads to an aperiodic sequence of propa- 
gation constants: e a et e a e a £b £a e&- • ■, where the first 
value corresponds to the value at the boundary of our 
semi-infinite array. We can redefine e a and tb using an 
appropriate adimensionalization in term of a single pa- 
rameter e (see Sec. II. C). In practice, since our array is 
of finite length N >■ 1, we only need to apply the substi- 
tutional rule a finite number of times. In the absence of 



B. Incommensurately modulated arrays 

Another way to generate a linear quasiperiodic array is 
by means of a spatial modulation of the refractive index 
that is incommensurate with the natural period of the 
array. An interesting example of this is the Aubry-Andre 
( AA) model [7 , where the propagation constant of each 
waveguide is given by 

£n = eo + V cos(27rnx) , (2) 



FIG. 2. (Color online). Intensity distribution of the linear 
modes for (a) periodic array, i.e., e = in both models, (b) 
Fibonacci array for e = 1, (c) and (d) refer to the A A array 
before (e = 1) and after (e = 3) the localization transition, 
respectively. The blue circle marks linear surface states. Only 
the half of each mode is shown to aid the visualization of the 
surface modes. Parameter v labels the eigenstates in increas- 
ing order of eigenvalue (mode propagation constant). 

where eo is the background value, r\ is the strength of 
the modulation, and \ is the ratio between the period of 
the modulation and the natural lattice period. Thus, the 
array is incommensurate when \ is an irrational number. 
The linear (7 = 0) AA model displays a sharp localization 
transition at \rj/V\ = 2: For \rj/V\ < 2, all eigenstates are 
extended, while for \r]/V\ > 2, they are all localized (see 
Fig. 2). Moreover, in the localized regime, the eigenvalue 
spectrum is composed by isolated points forming a fractal 
structure such that its complement is a dense set [18] . 

C. Adimensionalization 

For simplicity, we render Eq. ([!]) dimensionless by defin- 
ing E n {z) = y/VfiAnffifM*, with £ = Vz. Equation Q 
reads now 

i-^-A n + (A n+l +A n - l )+5 n A n +sgn(-y)\A n \ 2 A n = Q (3) 

where S n = (e n /V) — <j) and positive (negative) 7 denotes 
the focussing (defocussing) nonlinearity case. By a judi- 
cious choice of <j>, one can describe both, The Fibonacci 
model and the AA model. For the Fibonnaci model we 
define <j> as e a /V obtaining S n = or e and, for AA model 
we define <fi — £q/V obtaining S n = ecos(27rnx), with 
e = (ei, — e a )/V and e = rj/V, as the relevant control pa- 
rameters for the Fibonacci and AA models, respectively. 
The value of these parameters measure the strength of 
the quasiperiodic potential. 

The power content P = J2 n I^w| 2 j * s a conserved quan- 
tity in both models, and useful to characterize the differ- 
ent families of nonlinear modes. 



III. SURFACE SOLITONS 

We look for stationary localized surface mode solutions 
of pi). We pose the ansatz A n (z) = A n e lXz , obtaining the 
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FIG. 3. (Color online). Power vs propagation constant 
for different families of nonlinear surface modes for both fo- 
cussing and defocussing nonlinearities, and e > 0. Curves 
with positive (negative) A correspond to the focussing (defo- 
cussing) regime. Solid (dashed) curves denote stable (unsta- 
ble) branches. Top: Fibonacci array for e = 3. Bottom: AA 
array for e = 1 (phase of linear extended states). 



coupled system of algebraic equations 

-AA„+(A n+1 +yl n _i)+<5 n A„+sgn(7)| J 4 n | 2 A„ = 0, (4) 

where A n is real. Equation Q possesses the use- 
ful staggered-unstaggered symmetry (A, 7, e, A n ) — > 
(—A, —7, — e, (— l) n A n ). This means that only a half of 
the A — e space needs to be explored. The regions in pa- 
rameter space where the nonlinear modes are equivalent 
are: {7 > 0, e > 0} U {7 < 0, e < 0} which we denote 
as Type I nonlinear modes, and {7 > 0, e < 0} U {7 < 
0, e > 0}, denoted as Type II nonlinear modes. For each 
type, regions with 7 > and A > are associated to un- 
staggered modes, while regions with 7 < and A < are 
associated to staggered modes. In the remaining of this 
section, we will fix e > and both signs of 7 to examine 
the nonlinear surface modes present in the system. Equa- 
tions Q are solved by a straightforward implementation 
of a multi-dimensional Newton-raphson method |10j . us- 
ing the uncoupled (anticontinuum) limit to start the it- 
eration process . The linear stability of each mode is 
examined in the usual manner, by introducing a per- 
turbed solution of the kind A„(£) = (A n + ip n (Q) e iAS , 
with \ip n \ <C 1, where ip n is the (complex) perturbation. 
After replacing in Eq.Q and after linearizing in ?/>„, we 
obtain the evolution equation for the linear perturbation: 

= ^ ~ V'n-(V'n+l+V'n-l)-Sgn(7)A 2 l (2V'„+0- 

Next, the eigenvalue spectrum for perturbations of the 
type ipniO = ^nexp^f^) is examined, and the stabil- 
ity (instability) of the mode is inferred from the absence 
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FIG. 4. Stable nonlinear surface modes for values of power 
and propagation constants marked in Fig(3] Left: Fibonacci 
model. Right: AA model. 

(presence) of complex values of the perturbation eigen- 
values {Q,}. 

Figure [3] shows power content P vs propagation con- 
stant A curves for different families of nonlinear surface 
modes that inhabit the different gaps of the Fibonacci 
and A A models, as well as their stability, for both the 
focussing and defocussing regime. We observe that these 
modes originate from two main mechanisms: (i) Tangent 
bifurcations, usually observed in the external gaps, and 
giving rise to one stable and one unstable mode. Here, 
the local environment is predominant and the actual size 
of the array does not play a significant role, (ii) There 
are also nonlinear surface modes that bifurcate from lin- 
ear surface modes (Fig. la), or from extended asymmetric 
modes (Fig. lb). The most interesting feature of Fig (3] is 
the strong asymmetry between the focussing and defo- 
cussing cases, for the same distribution {e n } distribution. 

In Fig. [4] we show some examples of surface modes for 
parameter values in different regions and marked in Fig|3] 
They look qualitatively similar in both, the Fibonacci 
and AA models. Modes resembling these have also been 
observed in semi-infinite binary nonlinear waveguide ar- 
rays [T9] . 



A. Total Internal Reflection (TIR) modes 

TIR modes are highly-localized surface modes whose 
propagation constant resides in the outermost external 
gap, where total internal reflection mechanisms domi- 
nate. Usually, they need a threshold power to be excited. 
We observe that, for Type I modes, this minimum power 
increases with an increase in |e| (Fig{5]) in both models. 
For Type II modes, however, the behavior is different for 
each model. While for the Fibonacci model, the power 
threshold decreases monotonically with an increase in |e|, 
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FIG. 5. Power threshold for excitation of TIR modes as a 
function of e > for focussing and defocussing regimes. Solid 
(dashed) curves correspond to the Fibonacci (Aubry-Andre) 
models. 7 > and 7 < are associated to Type I and II 
nonlinear modes, respectively. 



for the AA model, the threshold decreases initially, reach- 
ing a minimum at |e| rs 2.59, where the minimum thresh- 
old is Pth ~ 1-25, and increasing monotonically with |e| 
afterwards (Figj5|. 

A linear stability analysis reveals that all these TIR 
modes undergo a stability change when dP/dX = 0, 
where the bifurcation occurs. The stability of the modes 
seems also to associated with a different position of its 
'center of mass': In stable modes (Fig.3a,b,d,e), the cen- 
ter of mass is closer to the surface, while unstable modes 
have their center of mass about one layer below the sur- 
face. 



B. Bragg Reflection(BR) modes 

These modes reside in deeper internal gaps in-between 
mini-bands, where Bragg reflection effects predominate. 
In these regions, localization is the combined result of 
nonlinearity plus the scattering produced by the aperi- 
odic structure. These modes are less localized than the 
TIR modes (see Fig. 3c, f), and they do not necessarily 
need a power threshold to excite them. Modes with 
no threshold bifurcate from two types of linear modes: 
Those localized at the surface and associated to isolated 
points of the linear spectrum (Figla), or those extended 
asymmetric modes where its envelope decreases away 
from the surface (Figlb). In any case, the presence of 
nonlinearity shifts the propagation constant of the linear 
mode outside the mini-band, giving rise to a nonlinear 
surface mode. 

In the linear Fibonacci case, and for a relatively large 
number of sites there are a finite number of linear modes 
localized at the surface. This number is a function of 
the total number of sites, N. In particular, for N = 89 
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FIG. 6. (Color online). Output power fraction at initial 
site vs quasiperiodic strength after a propagation distance of 
^max = 20. First (second) row: results for Fibonacci (Aubry- 
Andre) arrays and, left (right) columns: results for edge 
(bulk) excitation. Light and dark regions denote high and 
low degree of localization. The white curve denotes the power 
threshold to excite a TIR stationary surface mode shown in 
Fig|5] (only for surface cases). The dynamics of the marked 
states is shown in Fig[7] In bulk cases we averaged over 3 
realizations with different initial excite sites. 



there is a single localized linear surface mode (Figla), 
which gives rise to the families of Type I and II shown in 
Fig. 2. On the other hand, for the Aubry-Andre case, all 
of the linear eigenfunctions are extended for |e| < 2 and 
localized for |e| > 2[B], and some of them are asymmetric 
with maxima at or near the surface(Fig.lb). After the 
localization transition, these modes do not necessarily be- 
come localized at the surface (Fig.lc), but there are other 
linear modes that will become localized at the surface. 
When nonlinearity is present, we note that solitons re- 
siding in the internal gaps, between internal mini-bands, 
suffer a kind of oscillatory instability when their propa- 
gation constant approaches the edge of a mini-band, and 
their power content increases. This is due to a resonance 
between a linear mode at the edge of the band and the 
tail of the nonlinear mode, leading to a derealization of 
the nonlinear mode and to a synchronization of its phase 
with the one of the linear mode. Similar behavior has 
been reported in continuous systems with periodic po- 
tentials [2U] , and also in discrete binary systems [THl HI] ■ 



IV. DYNAMICS 

We look now into the effects of quasiperiodicity on the 
dynamical evolution of an initially localized input beam 
at the edge of the array and also inside the bulk, for 
comparison purposes. We solve Eq.Q numerically with 
an initial condition of A n (0) = \f~P S no . n , with uq = 1 for 
the edge excitation. For the bulk case, we have different 
local 'environments' around the initial site uq depending 
on the specific chosen site, we average the selftrapped 
fraction obtained for three different initial sites. 

For a given e and P, the space-averaged fraction of 
power remaining at the initial waveguide /, after a lon- 
gitudinal propagation distance £, m ax is defined as / = 

(P Cmax)' 1 Jq™° |-<4i(£)| 2 ^£- Results for / as a function 
of e and P with 7 > for both models, are shown in 
Figj6]in the form of a density plot. 

The first thing we notice from FigjH] is that both, Fi- 
bonacci and AA models, show a qualitative similar self- 
trapped fraction as a function of the aperiodic strength. 
The most interesting feature, however, is that the thresh- 
old power to excite localized modes is actually lower at 
the surface than at the bulk. More specifically, we note 
that for a given power content, a smaller quasiperiodic 
strength is required to effect localization at the surface 
than in the bulk. And for fixed quasiperiodic strength, a 
smaller power is needed to localize the excitation at the 
edge than inside the bulk. This is in marked contrast 
with the behavior observed in periodic systems. 

Coming back to the surface excitation case, we note a 
correspondence between the threshold power to excite a 
discrete surface soliton dynamically, with the thresholds 
for the stationary modes (white line) of the previous sec- 
tion. The dynamical thresholds are always larger. It is 
also interesting the marked asymmetry between the re- 
gion to the right of e = (Type I nonlinear modes) with 
the region to the left of e (equivalent to Type II non- 
linear modes). This asymmetry mirrors the one found 
for stationary surface modes. As expected, in the region 
where nonlinearity dominates, the output power fraction 
at the surface guide is close to unity, while at low input 
powers and high quasiperiodic strength, there is a sort of 
weak localization near the surface. At very small input 
power, the beam tends to diffract below certain quasiperi- 
odic strength. Conversely, when quasiperiodicity is very 
weak, the beam diffracts below a certain minimum in- 
put power. Different localization regimes are separated 
by regions where it is not possible to localize the wave 
function. Eventually, when the power or the strength of 
the quasiperiodic potential (e) is sufficiently high, most 
of the optical power remains trapped. 

Some examples of the detailed dynamical evolution in 
different regions of Fig[6]for surface excitations, is shown 
in FigjT] We see different types of localized surface solu- 
tions: 

(a) High-power solutions that are strongly localized at 
the surface and are associated to the TIR gap (Fig. 6a 
and 6e). Here, localization increases with power. 
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FIG. 7. (Color online). Dynamical evolution of a delta-like initial input beam, in the form of a density plot showing the 
amplitude |A„(;z)| 2 as a function of longitudinal distance z and waveguide position n. Plots (a)-(d) and (e)-(f) correspond to 
the Fibonacci and AA models, respectively. {z max — 20, N = 89). 



(b) Surface solutions that oscillate persistently between 
two sites, with a propagation constant that depends on 
the value of e (Fig. 6b and 6h). We have computed the 
Fourier transform of these solutions and found that the 
propagation constant associated to the transversal dy- 
namics (secondary propagation constant) is smaller than 
the propagation constant associated to the longitudinal 
dynamic (primary propagation constant), which domi- 
nates the dynamics. An increase in |e| render these type 
of solutions more localized and the secondary propaga- 
tion constant becomes negligible. However, their power 
content is smaller than the power associated with the 
TIR modes. This type of solution is similar to the ones 
described in Ref. [22) for a lattice with alternating site 
energies, where stable dynamic exists with the power os- 
cillating between two sites. 

(c) Low-power solitons that show oscillations with two 
different propagation constants (Figs. 6c, 6d and 6g): 
The first oscillation is fast and occurs between two nearby 
sites (as in case (b)); the second oscillation is slower 
and involves sites farther away from the surface (about 5 
sites). 

(d) Chaotic oscillations (Fig. 6f), where the oscillations 
do not seem to follow a well-defined pattern. These 
solitons reside at the very boundary between weak and 
strong selftrapped output fraction in the Pq — e diagram 
in Fig© 



V. NONDIAGONAL QUASIPERIODICITY 

For completeness, we also report here the results ob- 
tained for a variant of the Fibonacci model, where 
the aperiodic sequence occurs for the coupling between 
nearest-neighbor guides, and the individual refractive in- 



dex contrast of the guides are taken as identical. In that 
case, it can be absorbed as a phase factor in the solution, 
which we take equal to zero, without loss of generality. 
The equations of motion are: 



i~A n + V n+1 A n+1 + V n A„_ 1 + -/\A n \ 2 A n = 



(6) 



where V n denotes the coupling between guides n — 1 and 
n. According to the Fibonacci substitutional rule, The 
{V n } sequence is V a V b V a V a V b V a V b . . . 

The results obtained (not shown here) are qualitatively 
similar to the previous case of quasipcriodicity in the 
refractive index (diagonal quasiperiodicty) , with a lin- 
ear spectrum whose structure is composed of bands and 
gaps, with each band in turn, decomposed into smaller 
mini-bands and gaps. As in the previous case, nonlin- 
ear surface modes of the TIR and BR type can also 
be excited. For this diagonal quasiperiodic case how- 
ever, the usual staggered-unstaggered symmetry holds: 
(A, 7, A n ) -> (-A, -7, (-l) n A n ) as in the usual DNLS 
equation. This implies that the nonlinear surface modes 
for the focussing and defocussing cases are now equiva- 
lent. 



VI. DISCUSSION 

We have examined the interplay among nonlinear- 
ity, quasiperiodicity and surface effects in a semi- 
infinite, one-dimensional array of weakly-coupled, non- 
linear (Kerr) optical waveguides. We have considered 
two implementations of quasiperiodicity, namely, the Fi- 
bonacci and the Aubry-Andre model, and have found 
different families of surface localized modes, as well as 
the dynamical evolution of an initial surface input beam. 
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Both models show qualitatively similar behavior in the 
nonlinear regime, characterized by high-power solitons 
(TIR gap) strongly localized at the edge of the array, and 
lower-power solitons (BR gap) that oscillate persistently 
between some few sites close to the edge. Perhaps the 
main result is the strong asymmetry observed between 
the focussing and defocussing cases for the same distri- 
bution {e n }, that could be measured experimentally us- 
ing techniques similar to the ones used in ref [5]. Thus, 
for a nonlinear focussing medium, the defocussing case 
could also be explored by a judicious choice of the ape- 
riodic strength parameter e. For instance, in the Aubry- 
Andre model, where the modulation of the lattice is given 
by 7ycos(27m£), if the choice {77 > 0,7 > 0} gives rise 
to focussing modes, then to explore defocussing ones, 
is enough to change r\ — > —77, which in turn, is equiv- 
alent to explore the modes in the defocussing regimen 
{77 > 0,7 < 0}. A comparison of the dynamical evolu- 
tion of excitations initially localized at the surface and at 
the bulk, revealed that for a fixed degree of quasiperiodic 



strength, a smaller amount of optical power is required 
to selftrap the excitation at the surface than at inside 
the bulk. Conversely, for a fixed optical power, a larger 
quasiperiodic strength is required to selftrap the excita- 
tion inside the bulk than at the surface. To summarize, it 
is easier to localize an excitation at the edge than inside 
the bulk. This is in marked contrast with the behavior 
observed for periodic photonic lattices, and cast a note of 
caution on the widespread notion that a one-dimensional 
surface acts in a 'repulsive' manner. 
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